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A numerical study of the r-mode instability of rapidly 
rotating nascent neutron stars 
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ABSTRACT 

First results of numerical analysis of classical r-modes of rapidly rotating com- 
pressible stellar models are reported. The full set of linear perturbation equations 
of rotating stars in Newtonian gravity are numerically solved without the slow 
rotation approximation. A critical curve of gravitational wave emission induced 
instability which restricts the rotational frequencies of hot young neutron stars is 
obtained. Taking the standard cooling mechanisms of neutron stars into account, 
we also show the 'evolutionary curves' along which neutron stars are supposed 
to evolve as cooling and spinning-down proceed. Rotational frequencies of 1.4Af Q 
stars suffering from this instability decrease to around 100Hz when the standard 
cooling mechanism of neutron stars is employed. This result confirms the results 
of other authors who adopted the slow rotation approximation. 
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1 INTRODUCTION 

Recently the r-mode instabilit y induced by gr a vitational 
wave emission w as discovered (Andersson 1998; Friedman 
&: Morsink 1990- It is an axial mode variant of the secular 



instability of st ellar oscillations excited by coupling with a 
radia tion field ( Chandrasekhar 1970; Friedman & Schutz 



1978) 



The r-mode oscillations are the extension of Rossby- 
Haurwitz wave of a rotating fluid known in geophysi- 



cal literature s to global stellar oscillations (Papaloizou 



Pringle 1978). As they are dominated by the tangential 
velocity perturbations, their coupling with gravitational 
radiation was once considered to be weaker than those 
of the polar modes (e.g. f-modes). Thus it was so hI lik- 
ing and stimulating to the astrophysir-al i-ommunity that. 



enough to limit severely the rotational frequencies of hot 
young neutron stars which are born with the initial tem- 
perature of ~ 10 11 K. Even if the neutron stars are born 
with the nearly Keplerian (mass-shedding limit) rotation 
frequency, they settle down to states with only a small 
fraction of their original angular velocity, i.e. with less 
than 10 % of the Keplerian rotational frequency, when 
they cool down to the neutron superfluid transition tem- 
perature, T c ~ 10 9 K. As for the relatively cold and old 
neutron stars as those in X-ray binaries or radio pulsars, 
this instability ma y also set a severe limit on their rota- 

Andersson et al. 1999b|; 



tional frequ encies (Bildsten 1998; 



Levin 199£ ) 



Thus far, these results have been obtained by extrap- 



olating the results of the sl ow rotation approximation ( Pa- 



ma jiv people began to study whether or not the insta- paloizou fc Pringle 1978; ; |Provost et al. 198l| ; |Saio 198$ 



bility may be strong enough to set a severe limit on the 
rotationa l period of neutron stars. In particular - several 
authors ([Lindblom et al. 1998t Pwcn et al. 1998; Andcrs- 



al. 1999 



jon et al. 1999a; |Lockitch fc Friedman 1999; |Lindblom et 



to rotating models with the Keplerian frequency, except 
that Lindblom and Ipser (1999) obtained analytically the 
eigenfrequencies and the eigenfunctions of the 'classical' 
as well as the 'generalized' r-modes for the Maclaurin 



Yoshida fc Lee 2000] ) have investigated its effect 
on the evolution of the stella r rotational frequencies (see 
Friedman & Lockitch (2000) for a recent review). Those 



spheroids. Moreover, Bildsten & Ushomirsky (2000) re 



works indeed have revealed that the instability is strong 
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cently studied a damping effect of a viscous boundary 
layer formed at the stellar crust-core interface. Accord- 
ing to their work, the viscous boundary layer stabilizes 
the system significantly and the upper limit of the ro- 
tational frequency of neutron stars may become around 
~ 500Hz. If this is the case, the mode analysis without 
the slow rotation approximation should be needed. 
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Thus in order to see whether the evolutionary picture 
mentioned above is quantitatively correct and to investi- 
gate the evolution of stellar spins, it is indispensable to 
study the r- modes by taking rapid rotation into account. 

In this Letter, we will show the first numerical re- 
sults of the 'classical' r-modes for rapidly rotating stars. 
Here classical r-modes are oscillations with I = m, where 
I and m are the indices of spherical harmonics. Our nu- 
merical scheme is the imp roved version of the one us ed in 
the polar mode analysis (Yoshida & Eriguchi 1995) and 



will be published elsewhere together with mor e detailed 
res ults f or rapidly rotating compressible stars (Karino et 
al. Eoool ). 



2 EVALUATION OF DAMPING TIMES 

Once eigenfrequencies and eigenfunctions are obtained, 
timescales of the change of the system due to gravitational 
radiation a nd viscosity can be es timated by the standard 
procedure (Lindblom et al. 1998). 



The canonical energy E c of the perturbation is de- 
fined as: 



d x ■ [pSv ■ Sv * + (8p/p — 8(f>)*8p\ 



(1) 



where 5p, Sp, Sv, 5<f> are the Eulerian perturbations of the 
density, pressure, velocity and gravitational potential, re- 
spectively. The superscript, *, means the complex conju- 
gation of the corresponding quantity. 

When we assume that the perturbed quantities be- 
have as ~ e -*(wt-n»¥>) ^ thg energy dissipation rate due to 
gravitational radiation is expressed by the multipole ra- 
diation formula as follows: 



E gr = (id — mft) y NlLO 



l>m 



2Z + 1 



(\6D, 



/„.| 2 + \SJhn\ 2 ) 



(2) 



where SDi m and SJi m are the mass and the mass current 
multipoles, respectively, defined as: 



SD tm - J d 3 x- Spr l Y* m , 



SJi, 



l + l 



d x ■ r (pSv + 5pv) ■ Y lr 



(3) 



(4) 



where c is the speed of light. Here the magnetic type vec- 
tor spherical harmonics Y^ is defined as: 

?& = [i(i + i)]-$i>xvy, m . (5) 

The coupling constant Ni is a function of I as follows: 



Nt 



4wG (Z + l)(Z + 2) 
c 2 '+U(Z-l)[(2Z + l)!!] 2 ' 



(6) 



As the rapidly rotating stars are highly deformed 
from spherical configuration, the summation (2) in the in- 
dex I should include theoretically the infinite numbers of 
multipole terms. We take the lowest five ones into account 
here for simplicity. It should be noted that the slow rota- 
tion approximation includes only the lowest order term. 
On the other hand, the energy dissipation rate due to the 
fluid viscosity is written as: 



E s = - / d s x ■ 2riSa ab 5a* a 



(7) 



for the shear viscosity with a coefficient r\, and 

E h = - I d s x-C\SQ\ 2 , (8) 

for the bulk viscosity with a coefficient (. Here the volume 
expansion rate SO and the shear tensor 8a a b are defined 

as: 

(9) 



se = v c <5v c , 




and 




8a ab = V( a <5u b ) - 


se 

- -17-Oab, 



(10) 



where a round bracket means symmetrization among in- 
dices. These energy dissipations and their rates are nu- 
merically evaluated for the rotating stellar models. 

Then the imaginary part of the eigenfrequency of the 
r-mode, t" 1 , is written as follows: 



-i _ E gr + E s + Et _ _i _i _i 

T r — ^ = T gr + T s + T b ■ 



(11) 



Positive values of r" 1 correspond to the states where the 
r-mode instability dominates the stabilizing viscous effect. 
Consequently critical states of the instability are defined 

by t-,7 1 = 0. 



3 COMPARISON WITH THE 

SLOW-ROTATION APPROXIMATION 

We have studied the adiabatic perturbations of the equi- 
librium sequence of 'canonical neutron stars', i.e. poly- 
tropes with p — Kp 2 whose mass and radius in the spher- 
ical limit become IAMq and 12.5km, respectively. Here 
p, p and K are the pressure, the density and a polytropic 
constant, respectively. This sequence is the same as that 
analyzed by many authors in the slow rotation approxi- 
mation. In Fig 1, the eigenfrequency u) normalized by the 
angular velocity £1 is plotted against the normalized an- 
gular velocity fi/ ' ^JAiiGpc where p c and G are the central 
density and the gravitational constant, respectively. f\ To 
compare our results with those by the slow rotation ap- 
proxi mation, we also plot the solutions by Yoshida and 
Lee (Yoshida fc Lee 200C|) where the third order term 



in the f2 expansion is taken into account. The result of 
the slowest rotation in our results (£1/ \f4-n Gp c — 0.0154) 
agrees with theirs within the relative error of 0.4%. The 
grid point numbers used in this computation are 32 for the 
radial coordinate r from the center to the surface of the 
star, and 11 for the angular coordinate 9 in the quarter 
of the meridional cross section of the star. 



T We do not display the results of our computation in the vicin- 
ity of the mass-shedding limit, because the eigenfunctions are 
rather ill-behaved and the estimate of the instability time-scale 
is difficult there. Our results (Fig. 1 and Fig. 2) are terminated 
at the model with 98.7% of the mass-shedding limit frequency 
(850Hz). 
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n/(4wGp c ) ' 



Figure 1. The normalized eigenfrequcncy tu/Q is plotted 
against the normalized angular velocity r2/(47rGp c ) 1 ' 2 . Also 
plotted are the solutions by the slow rotation approximation. 
The dashed line is the solution of the first order in f!, whereas 
the dotted line is the one in which the third order correction 
is also taken into account. The mass-shedding limit of this se- 
quence is the point where f2/(47rGp c ) 1 ' 2 = 0.205 




Temperature (K) 



Figure 2. Critical rotational frequency of a IAMq neutron 
star with the polytropic equation of state p = Kp 2 is plotted 
against the temperature by the thick solid curve for the I = 
m = 2 mode. For comparison the critical curve computed using 
the slow rotation approximation is drawn as well (thin solid 
line). Also shown are the curves on which the cooling time- 
scale of the star is equal to the r-mode instability time-scale 
(dotted line for the modified URCA process; dashed line for 
the 'quasi-particle' /3-decay in pion condensation). 



4 CRITICAL CURVE AND 

'EVOLUTIONARY TRACK' OF A HOT 
AND YOUNG NEUTRON STAR 

As for the source of shear and bulk viscosity, we adopt 
the same microphysical process as that used by Ipser and 
Lindblom ([pser fc Lindblom 1991); that is, if the stellar 



temperature T satisfies the relation T > T c , the dominant 
processes are neutron-neutron (nn) collision and electron- 
electron (ee) collision. They result in the shear viscosity 
coefficient with the following density and temperature de- 
pendence: 



9/4 T - 



Vnn OC p 

and 



7? ee oc p 2 T 2 . 



(12) 



(13) 



When T < T c , only the electron collision is included. The 
bulk viscosity arises from lag of the /3-reaction to the os- 
cillation. The viscosity coefficient has the following de- 
pendence: 



£ oc p (lu — mQ.) 



-2m6 



(14) 



For each mode, the inverse of the time-scale r" 1 de- 
pends on the stellar temperature T as well as on the stel- 
lar rotational frequency /. In Fig El we plot the critical 
points of stability on the T — f plane for the I = m = 2 
mode. Above this critical curve, gravitational wave emis- 
sion instability dominates the stabilization effect due to 
viscosity. Along the fixed temperature line, i.e. vertical 
line, the instability time-scale decreases approximately as 
~ y 2m + 2 . As seen from this figure, the smallest value of 
the rotational frequency of the critical curve is roughly 8 
% of the frequency at the mass-shedding limit. 

Also shown is the critical curve using the slow rota- 
tion approximation. The expression of the timescales is 
adopted from Lindblom et al. (1999). It is noted that the 
shear viscosity term (Eq.(10)) is evaluated in our analysis 
by numerical differentiation of the corresponding eigen- 
function. This leads to the relative error of ~ 10 % in the 
shear viscous timescale r s for the present mesh size (32 
and 11 in the r and 9 direction). 

Together with it plotted are the curves on which in- 
stability time-scale is equal to the cooling time-scale of the 
star. If we assume the 'stan dard' modified URCA process 
( Shapiro fc Teukolskyl983| ) to be dominant in the initial 
stage of neutron star cooling^] the effective cooling time- 
scale t coo i at the temperature T(K) is defined by: 



7"cool 



rilllTg 



dt 



6tc 

1 9 



(15) 



where Tg = T/10 and t c is a constant characterizing the 
cooling time which is typically ~ ly for this process. 

Apart from the early phase, hot young neutron stars 
may evolve along this evolutionary curve r r = T" CO oi- Sup- 
pose that a neutron star is born at the temperature 
~ 10 11 K and with the rotational frequency above the 
bottom of the critical curve. At the beginning it is sta- 
ble against the r-mode instability and evolves along the 
horizontal line on the T — f plane. In a few second, it 
enters the unstable region against the r-mode and the ini- 
tial (stochastic) perturbation begins to grow. As long as 
t t > Tcool, it evolves almost horizontally. Once it reaches 

-f Though nucleon bremsstrahlung is also one of the stan- 
dard cooling channels, its contribution is order of magnitude 
smaller tha n that of the modified URCA process. See Friman 
& Maxwell (Il979|). 
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the 'evolution curve', the star begins to go down along 
this curve, since any hypothetical displacement of the star 
from this curve will be amended by the growth of insta- 
bility or the cooling. If the star is located at the right 
side of this curve, the star will evolve nearly horizontally 
leftwards due to rapid cooling. If the star is located at the 
left side of this curve, the star will evolve nearly vertically 
downward due to rapid loss of gravitational waves. 

When the star goes down along the evolution curve 
and approaches the critical curve asymptotically, the in- 
stability ceases to work effectively and the star will not 
be spun down by the instability any more. The location 
of this evolution curve in the T — f plane is insensitive 
to the value of constant t c . Consequently every neutron 
star seems to settle down to an almost universal terminal 
rotational frequency, i.e. ~ 100Hz at the end of the spin 
evolution. If t c — ly, this frequency is 90Hz and is reached 
when the temperature of the star equals to 1.1 x 10 9 K. 

On the other hand, if some exotic rapid cooling mech- 
anisms dominate, the picture may be changed drastically. 
In Fig 0, an evolution curve defined by r r = r coo i for 
the model with pion condensation is also plotted (dashed 
curve). The pion condensation enhances the cooling be- 
cause an^w_Drocess_ariaj£g_ojis_to the ordinary URCA pro- 
cess ( phapiro fc Teukolskyl983 ) begins to operate. The 
effective cooling time-scale for this process is written as 
follows: 



Tcool 






(16) 



where t' c ~ 200s. The evolution curve of r r = r coo i shifts 
upward and the terminal rotational frequency also shifts 
to about 200Hz. 

Commenting on an evolution of a star by a rapid 
cooling mechanism may be appropriate here. By an ex- 
otic cooling mechanism such as the pionic reaction, the 
core fluid of a hot and young neutron star cools rapidly 
down to ~ 10 8 K, while the crust component of it remains 
hot (~ 10 9 K) in the first few decades (see Lattimer et 
al. (1994) for instance). Then the standard assumption 
in the r-mode instability that the stellar te mpera ture re- 
mains uniform is broken. Andersson et al. (1999a) argues 
that the warm crust work as a heat reservoir and the 
rapid cooling effect on the stellar evolution by the r-mode 
is moderated. This reduces the difference between the out- 
come of the rapid and the standard cooling. We further 
note that th e viscous heating by the r-mode shear motion 
( |Levin 1998| ) , which is omitted here and in the other works 
dealing with hot and young neutron stars, may delay the 
cooling and have the same effect as the crust. This heat- 
ing prolongs the lapse during which the r-mode instability 
is effective (Yoshida, unpublished). 

This means that in order to know the significant ef- 
fect of the r-mode instability on spin evolution of newly 
born hot neutron stars, we need to know more about the 
nature of dense matter at high and intermediate temper- 
atures. 



5 SUMMARY 

We show the first results of the numerical mode analysis 
of the classical r-mode for the rapidly rotating and com- 
pressible stellar models. The eigenfrequencies and their 
eigenfunction of the rapidly rotating stellar models are 
numerically computed nearly to the Kepler limit. The 
time scale of the r-mode instability is evaluated by the 
eigenfrequency and the eigenfunction. 

The behaviour of the eigenfrequency is confirmed to 
be well approximated by the third order expansion in the 
angular frequency fi. The difference between the critical 
curves drawn by the full computation and by the slow 
rotation approximation is shown to be rather small. This 
certificates the use of the slow rotation approximation in 
the third order to investigate the r-mode instability. 
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